Quantum entanglement in the NMR implementation of the Deutsch-Jozsa algorithm 
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A scheme to execute an n-bit Deutsch-Jozsa (D-J) algo- 
rithm using n qubits has been implemented for up to three 
qubits on an NMR quantum computer. For the one and two 
bit Deutsch problem, the qubits do not get entangled, hence 
the NMR implementation is achieved without using spin-spin 
interactions. It is for the three bit case, that the manipula- 
tion of entangled states becomes essential. The interactions 
through scalar J-couplings in NMR spin systems have been 
exploited to implement entangling transformations required 
for the three bit D-J algorithm. 

The utilization of the intrinsically quantum mechanical 
nature of the physical world to widen the scope of com- 
putational algorithms is one of the important discoveries 
of this decade @-§. It was shown recently that quan- 
tum computers can perform certain computational tasks 
exponentially faster than classical computers 

The Cleve version of the D-J algorithm which re- 
quires n + 1 qubits to solve the n-bit Deutsch problem, 
has been implemented by several research groups using 
NMR js] [ll] . It has been shown recently that the n-bit 
Deutsch problem can be solved using n qubits alone [jl2| . 

In this paper, we experimentally demonstrate that the 
n-bit D-J algorithm does not require n + 1 qubits for its 
implementation. By doing away with the extra qubit, 
the algorithm can be more easily accessed for a greater 
number of qubits. Furthermore, the one-bit and two-bit 
implementations of the modified version do not involve 
quantum entanglement |I^,[l3|]. In these cases, only the 
concept of coherent superposition is exploited, to pre- 
pare "in parallel" an input state which is a superposition 
of all possible classical inputs, and the experiment has 
been performed without using spin-spin interactions. It 
is only in the implementation for three or more qubits 
does quantum entanglement play a vital role. A judicious 
combination of rf pulses and free evolution intervals (un- 
der the interaction Hamiltonian), has been employed to 
construct the required entangling transformations. 

Consider an n-bit binary string function / can 
be defined on this n-bit domain space to a 1-bit range 
space, with the restriction that either the output is the 
same for all inputs (the function is constant) or the out- 
put is for half the inputs and 1 for the other half (the 
function is balanced). All the 2 n possible input strings 
are valid inputs for the function (f(x) = {0, 1}). In quan- 
tum computation, these n-bit logical strings are in one- 
to-one correspondence with the eigen states of n-qubits, 
and one can hence label the logical string x by the eigen- 
state | a;). Classically, for an n-bit domain space, one 
needs to compute the function at least 2™/ 2 + 1 times in 



order to determine whether it is constant or balanced. 
The D-J algorithm achieves this on a quantum computer 
using only a single function call [^|,[7j . 

The usual implementation of the D-J algorithm for n 
bits requires n + 1 qubits, the function / being encoded 
through an /-dependent unitary transformation, 
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where © denotes addition modulo 2. The implementa- 
tion of the unitary transformation Uf, along with the 
Hadamard transformation, then suffices to distinguish 
the function as constant or balanced j|]7|]. A Hadamard 
transformation on one qubit mixes the eigenstates max- 
imally, 
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The Hadamard transformation for n-qubits is the tensor 
product of the one-qubit transformation (H n = H (g> H £g> 
) H), its action on the n-qubit eigen states being 
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where Xj and yj are the jth entries of the n-bit strings x 
and y. 

A modified scheme can be designed to solve the li- 
bit Deutsch problem, using n qubits alone p"2j . Here, for 
every function / a unitary transformation is constructed, 
such that its action on the eigenstates of n-qubits is 

|*U t (-l) f ^\xU lt (4) 

Consider n qubits, all in the state |0) ; a Hadamard 
transformation H n converts this state to a linear su- 
perposition of all 2™ eigenstates with equal amplitudes 
and no phase differences. The unitary transformation 
Uf (defined in Eqn. ||) acting on this state, introduces 
an /-dependent phase factor in each eigenstate in the 
superposition. At this juncture, all information about 
/ is encoded in the quantum state of the n qubits. A 
Hadamard transformation H n is once again applied in 
order to extract the function's constant or balanced na- 
ture: 
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The final expression for the output state in Eqn. || has 
an amplitude 1 for the state |0) n _ bit for a constant func- 
tion and an amplitude for a balanced function. The 
categorisation of the function as constant or balanced 
through a single function call using n qubits, is shown 
pictorially in Fig. [I]. 
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FIG. 1. The block diagram for the modified D-J algorithm. 

The number of functions for the n-bit Deutsch prob- 
lem is n Cn/2 + 2 (where N = 2"). The experimen- 
tal implementation of the modified D-J algorithm for n 
bits requires the realisation of the unitary transformation 
corresponding to each of these functions, and the n-bit 
Hadamard transformation, on a physical system. We now 
proceed towards the NMR implementation of the modi- 
fied D-J algorithm for one, two and three qubits, where 
the number of functions are 4, 8, and 72 respectively. 

The pseudo-Hadamard transformation [lj] (practi- 
cally equivalent to the Hadamard operator described in 
Eqn. ||) achieved by a (90° ) y pulse on a spin, has been 
utilised in our experiments. For the case of two and three 
qubits, the same has been achieved by a (90°) y pulse ap- 
plied non-selectively on all the spins. 

The n-bit unitary transformations Uf corresponding to 
the functions /, are diagonal in the eigenbasis and find a 
natural description in terms of the single-spin operators, 
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where j labels the qubit involved. The action of Uf on 
an eigenstate (as described in Eqn. can been used 
to calculate the explicit matrix forms of Uf, for every 
function /. 

The operator representations of the four unitary trans- 
formations for the one-bit modified D-J algorithm are 



(7) 



A pseudo-Hadamard operation achieved by a (90°)j, pulse 
is applied on a thermal initial state, in order to cre- 
ate a coherent superposition prior to applying the de- 
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sired unitary transformations Uf. The constant func- 
tions correspond to a "do-nothing" operation, while the 
balanced functions are achieved by a rotation by the an- 
gle 7r about the z-axis of the spin, upto a global phase 
factor. These z-rotations have been implemented using 
composite-z pulses, whereby a rotation by an arbitrary 
angle 9 about the z-axis, can be decomposed as a set of 
rotations about the x and y axes pM. 
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Global phase changes are not detectable in NMR and are 
hence ignored. The spectrum reflects the constant or bal- 
anced nature of the function implemented (Fig. ||). After 
the implementation of a balanced function, the qubit is in 
a state out-of-phase with the rest of the spectrum. The 
modified D-J algorithm for one qubit demonstrates the 
power inherent in even a single bit of quantum informa- 
tion. 
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FIG. 2. The modified D-J algorithm for one qubit imple- 
mented on 5-nitro-2-furaldehyde, the proton resonating at 9.2 
ppm being chosen as the single qubit. 

All the 8 unitary transformations corresponding to the 
functions for the two-qubit case are given in terms of J( J ) 
and Gz 



as 
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All these operators are direct products of single spin op- 
erators. They are thus incapable of generating entangled 
states and can be implemented by operations on indi- 
vidual spins. A pseudo-Hadamard transformation was 
performed on all the spins (initially in thermal equilib- 
rium) prior to the execution of the desired Ui transfor- 
mations. The two constant functions correspond to the 
"do-nothing" operation, experimentally. The NMR im- 
plementation of the balanced functions U^ ~ blt) and U^' bA) 
involves rotations by 7r about the z-axis in the single- 
spin subspaces of spins 1 and 2 respectively, and have 
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been achieved using composite-z pulses (Eqn. ||). The 
^r(2-b,t) transformation too, does not require the scalar J 
spin interaction and is implemented as successive tt rota- 
tions about the z-axes of spins 1 and 2 respectively. The 
balanced functions are distinguished by one (or both) the 
spins being out-of-phase with the rest of the NMR spec- 
trum (Fig. ||). Only half the total number of functions 
have been shown in the one and two-qubit cases, as the 
others are merely negatives of these, and lead to the same 
spectral patterns (spectra not shown). 
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FIG. 3. The modified D-J algorithm for two qubits, imple- 
mented on 5-nitro-2-furaldehyde with the two qubits resonat- 
ing at 6.47 ppm and 6.29 ppm respectively. 

The three-qubit D-J algorithm affords the simplest ex- 
ample where quantum entanglement plays a definitive 
role in the computation. The task here is to implement 
72 unitary transformations, the explicit operator forms 
for 9 of which are 
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The operators {/ ; J 3 " b, * ) - ^ 3 b,t ' can be decomposed as 
direct products of single-spin operators and are thus 
non-entangling transformations. The operators £/g 3 ~ blt> 



- Ug' blt} cannot be decomposed as direct products of 
single-spin operators and are hence capable of generat- 
ing entangled states from non-entangled ones. The op- 
erators Uq , Uf' blt) , and t/g 3 blt) , are entangling in dif- 
ferent two-spin subspaces and can be factored as direct 
products of a single-spin operator and a two-particle en- 
tangling transformation. They can thus generate states 
in which two qubits are entangled, with the third qubit 
remaining non-entangled with either of them. On the 
other hand, the transformation Ug 3 ' b,t> does not allow any 
such simplifications. It is maximally entangling and leads 
to states that are three-qubit entangled. These 9 func- 
tions are thus divided into three categories namely, non- 
entangling, two-qubit entangling, and maximally (three- 
qubit) entangling. The remaining functions are similar 
in form, and can be classed into one or the other of these 
categories. 

The result of experimentally applying the non- 
cntangling transformations JJ^' blt> - L7^ 3 blt) ; after a 
pseudo-Hadamard transformation on a thermal equilib- 
rium state, is shown in Figure ||. The constant func- 
tion £/.j 3 ~ blt) is the unit operator, and corresponds to the 
"no pulse" or the "do-nothing" operation. The balanced 
functions t/^ 3 blt) and t/3 3 ~ blt> correspond to a rotation by 
the angle tt about the z-axis of the first and the third 
spins respectively, without perturbing the other spins. 
This has been achieved by a spin-selective [tt] z pulse in 
each case, using composite z-pulses (Eqn. ||). The spec- 
trum is categorised by the spin in question being out-of- 
phase with the rest of the spectrum. The transforma- 
tion U^ > ' h ' t) has been implemented by two spin-selective 
[tt] z pulses applied consecutively on the first and the 
second spins respectively, and leads to a spectrum with 
both these spins being out-of-phase with the third. The 
non-entangling balanced function [/g 3 ~ blt) , has been im- 
plemented by successive spin-selective [tt] z pulses on all 
the three spins. 



The two-qubit entangling transformation Uq 



is 



achieved by the pulse sequence 

3 ( T23 / 2 ) ,5 



(1-23/2) 



[tt/2] 

where T23 = I/J23 and 1,2 and 3 are qubit labels. The 
operators Uf' blt) and U%~ blt) correspond to cyclic permu- 
tations of the qubits. The spin-selective tt pulses in the 
middle of the free evolution period T23 refocus the chem- 
ical shift evolution. The pulse sequence (applied after a 
pseudo-Hadamard transformation on all three qubits in a 
thermal initial state) results in a density matrix with the 
product operator form — 1\ + 2I 2 / 3 + 2J 2 / 3 , leading to 
a spectrum with the multiplet of the first qubit inverted, 
and an antiphase doublet of doublet pattern for the other 
two qubits (Fig. ^). The three-qubit entangling function 
[/g 3 ~ blt) is implemented by the pulse sequence 



[*/2]i [tt/2]] 
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leading to an antiphase spectral pattern for all three 
qubits, corresponding to the product operators 2J*J| + 
2/2/3 + 4/1/2/3 ( Fig ^ |) The spectra in Figure | suffer 

from phase distortions arising from the inaccurate refo- 
cusing of the chemical shifts during the r periods, and 
J-evolution during long spin-selective composite-z pulses 
(Ml = [^r/2]*,[7r]* [7r/2]*_ x = 42 msec). 
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well as a good initial state. The observable spectral re- 
sult is the same in both cases, though beginning with a 
pseudo-pure state creates some (undetectable) multiple- 
quantum coherences. The final pseudo-Hadamard trans- 
formation to extract the constant or balanced nature of 
the function (Fig. [I]) is canceled by the (90)° read-out 
pulse usually used in NMR experiments, and the compu- 
tation essentially culminates in the application of the de- 
sired Uf function after the first pseudo-Hadamard trans- 
formation. 

A modification to the usual D-J algorithm enabled 
an n-bit implementation using n qubits. The required 
unitary transformations were tailored to eliminate the 
need for the extra qubit, and the modified D-J algorithm 
was tested experimentally for one, two and three qubits. 
While the one and two qubit cases use non-entangling 
unitary transformations, it was noted that for three (or 
more) qubits, multi-particle entangling transformations 
are required. 
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FIG. 4. The refined D-J algorithm for three qubits, imple- 
mented on 2,3-dibromopropionic acid. The functions shown 
are all non-entangling in nature. 
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FIG. 5. Entangling balanced functions implemented on the 
three-qubit system of 2,3-dibromopropionic acid. 

This implementation of the D-J algorithm does not re- 
quire the initial preparation of the spins in a pseudo-pure 
state, since the thermal equilibrium state serves equally 
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